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0. Background 

This report is one of a series of studies steinming from 
current research on electrohydrodynamic (EHD) power genera- 
tion, sponsored by the Department of Energy. A certain EHD 
scheme proposed by the Marks Polarized Corporation is presently 
under development. The authors of this report have partici- 
pated as consultants on various theoretical and analytical 
aspects of this effort. One of our main objectives has been 
to make an independent appraisal of the performance possi- 
bilities and limitations inherent in the general type of EHD 
system involved in this development. Another objective has 
been to create a sound analytical basis for optimizing the key 
design parameters of such a device, and to suggest design 
improvements where appropriate. 

The results of our first study are reported in detail in 
Ref. (1) and are summarized more concisely in Ref. (2). This 
v7ork indicated that the design as initially proposed would 
yield unacceptably poor performance and suggested a certain 
design revision. Further work was done to analyze such a 
revised design and is reported in the appendix of Ref. (1) . 

An improved, simplified and shortened version of this work is 
reported in Ref. (3). Unfortunately, our results indicated 
that the revised design, while showing somewhat improved per- 
formance, still would not attain competitive overall efficiency. 

The Marks Polarized Corporation has challenged these 
pessimistic conclusions. They argue that in our analysis 
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the compressibility of the medium was not taken into account 
and that if it had been, the predicted performance would have 
been drastically improved. We have good reasons to conclude 
that Marks' argument about compressibility is fallacious. In 
the first place our analysis did consider all of the major 
effects of compressibility; only certain minor effects in one 
part of the system were neglected in order to simplify the 
calculations . While a more elaborate and accurate treatment 
of these neglected effects is certainly possible, we maintain 
that this refinement will not suffice to change the generally 
pessimistic conclusions established in our original analysis. 

Fortunately, this difference of opinion between ourselves 
and the Marks Polarized Corporation pertains to a question of 
fact about which it is entirely possible to obtain objective 
evidence. To this end we present in this report a revised 
and improved analysis which takes all compressibility effects 
exhaustively and minutely into account. Calculations based on 
this revised analysis should eventually settle the above 
question one way or the other. A typical trial calculation is 
presented later in this report. This initial calculation has 
been made manually, but systematic further work will be done by 
computer. The full results and conclusions derived from these 
calculations will be reported later. 

It has also been conjectured by Marks and others that it 
might be possible to improve greatly the performance of the 
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EHD generator by employing a two-fluid cycle. We therefore 
present in this report a detailed analysis and procedure by 
means of which the performance of any specified two-fluid EHD 
system can be calculated. Our analysis takes into account 
fully the compressibility of both fluids at all locations in 
the system. 

We propose to undertake systemic computer calculations 
of this kind which should enable an objective assessment to 
be made of the performance potentialities and limitations of 
several proposed two-fluid EHD generators. Such results would 
be of great value in helping decide whether the two-fluid EHD 
generator shows sufficient promise to warrant further develop- 
ment effort and, if so, what the optimum design point values 
of the key parameters should be. 

As far as we know, the present report is the only pub- 
lished analysis of the two- fluid EHD generator which includes 
the compressibility of both fluids and which considers the 
complete thermodynamic cycle. The principal earlier effort 
along these lines is given in Ref. (4) , but the present report 
goes well beyond anything undertaken in Ref . (4) . In particu- 
lar, Ref. (4) makes no attempt to analyze the complete thermo- 
dynamic cycle. Moreover, although it contains certain useful 
data and is correct in some respects. Ref. (4) is flawed by a 
number of serious errors . 
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While the present report is otherwise self-contained, it 
does assume that the reader has a general familiarity with the 
type of EHD power generator described in Refs. (1) , (2) and 
(3) . It is recommended that the reader unfamiliar with this 
general scheme first review briefly the introductory material 
in the first section of any one of these references before 
proceeding with the more detailed treatment in this report as 
this earlier background information is not covered or repeated 
in the present report. We do, however, include here a schema- 
tic diagrara. Fig. 0.1, which is taken from Ref. (3). It 
illustrates one particular version of an EHD power generator. 
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Introduction 



For purposes of analysis, the overall thermodynamic pro- 
cess in the two-fluid EHD generator may be represented by the 
schematic flow diagram shown in Fig. 1.1. 

Primary fluid enters the ejector at station 1, secondary 
fluid enters at station 2 and the mixture leaves at station 3. 

The mixture then passes through a diffuser, process 3-^4, 
the basic purpose of which is to reduce the kinetic energy of 
the flow and thereby decrease the associated friction losses 
downstream of the diffuser. On the other hand the diffuser 
introduces certain losses of its ov;n so that some care is 
necessary to optimize the design in this respect. 

The gas mixture then passes through the electrical power 
section, process 4->5 , from which the gross electrical power is 
extracted. 

Finally, the gas mixture enters the condenser/separator 
at station 5. The primary fluid is separated from the gas 
mixture by condensation and leaves at station 6. The fluid at 
this point is assumed to be compressed liquid at a known static 
pressure Pg and a known static temperature Tg . Kinetic 
energy at station 6 is regarded as negligible. The secondary 
fluid leaves the condenser/separator and enters the ejector at 
station 2. We assume, moreover, that 

Pg = Pj = p^ (1.1) 

and 

Tg = T 2 (1.2) 
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Fig. 1.1 SCHEMATIC FLOW DIAGRAM OF TWO - FLUID EHD CYCLE. 




The prima3ry fluid which leaves the condenser/separator 
at station 6 is circulated back through the pump and through 



the boiler/superheater, process 7^1. 

It should also be explained that the primary fluid at 
station 1 is a condensible vapor of high molecular weight 
which contains fine liquid droplets. Unless specifically 
stated otherwise the quality z at this location is taken 
as 0.95. The important thermodynamic static properties at 
station 1 can be expressed in the following form. 



Here subscript f denotes the saturated liquid and s\abscript 
fg denotes the change between saturated liquid and saturated 
vapor. All properties are evaluated at a specified pressure 
; the corresponding saturation temperature T^ is then 
also known. 

The corresponding stagnation properties T^^ and h^^ 
are fixed when the stagnation pressure is specified. 

Of course the entropy has the same value s^ at the stagna- 
tion condition as at the corresponding static state. Thus 
specification of P^^^ and s^^ suffices to fix all other 
stagnation properties. 

While it is convenient to start the analyses by stipula- 
ting the quality at station 1, it is also advisable later in 
the calculation to determine the corresponding quality at 



h^ = h^ + z h^g Enthalpy 



(1.3) 



s^ = s^ + z Sjg Entropy 



(1.4) 
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station 4, the entrance to the electrical working section. 

It is essential that the fluid entering the working section 
contain finely dispersed liquid droplets of the proper size. 
These droplets carry the electrical charges that are the 
central feature of EHD power generation. We assume tentatively 
that acceptable droplet characteristics are obtained by 
maintaining the quality at the inlet to the working section 
in the range 0.92 to 0.98. These limits are only estimates. 

The thermodynamic analysis would be greatly simplified 
if both fluids were perfect gases with constant specific 
heats. In fact only the secondary fluid satisfies this re- 
quirement. Nevertheless, we can define an "equivalent perfect 
gas" which adequately approximates the essential thermodynamic 
properties of the primary vapor/liquid mixture. When carried 
out judiciously, this procedure simplifies the analysis at the 

5 

cost of only a small loss in numerical accuracy . 

The ratio of specific heats y of the "equivalent perfect 
gas" may be deduced from the important relation 




Yi/(Yi - 1) 



(1.5) 



Solving this for gives 



^1 " £n(Pg^/P^) 



_ ^n(T^3^/T^) 



( 1 . 6 ) 



Another fundamental property of the "equivalent perfect 
gas" is its specific heat which may be evaluated from 

the expression 
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(1.7) 



pi " - T^) 

where all quantities on the right are now known. 

It then follows from standard perfect gas relations that 
the gas constant and the molecular weight of the 

"equivalent perfect gas" are, respectively, 



^1 = 






Si 



( 1 . 8 ) 



and 




where R = universal gas constant 
= 8315 joules/kg °K 

Eqs. (1.6) through (1.9) complete the definition of the 
"equivalent perfect gas" which adequately simulates the essen- 
tial thermodynamic properties of the primary vapor/liquid 
mixture . 

It should also be explained that the secondary fluid at 
station 2 is predominantly a noncondensing gas of low molecu- 
lar weight primary vapor. 

It can be shown that the mass ratio o of condensible 
primary vapor to dry secondary gas at station 2 is given by 
the expression 



"g IP2 - ‘ 



( 1 . 10 ) 
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where 



= true molecular weight of condensible vapor (not to 
be confused with the "equivalent molecular weight 
W^" , considered earlier) 

= molecular weight of noncondensing gas 
Pv(Tg) = vapor pressure of condensible primary fluid at known 
temperature Tg 

?2 = known static pressure at station 2 

The present analysis is restricted to conditions under 
which CO is very small compared with unity and may be neglec- 
ted. Thus the fluid at station 2 may be treated as dry gas. 
This greatly simplifies the analysis. Eq. (1.10) is useful for 
verifying that this assumption is indeed satisfied in any par- 
ticular instance. 

It should be added that if it later becomes of interest 
to make calculations for circumstances under which parameter 
CO is not negligible, the present analysis can in fact be 
generalized to include this effect. For the present, however, 
we prefer to deal only with the simpler situation in which co 
may be neglected. 
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2. Ejector 

The ejector receives a primary stream of high molecular 
weight gas at station 1, a secondary stream of low molecular 
weight gas at station 2 and discharges the resulting mixture 
at station 3. Static conditions at these three stations are 
designated by subscripts 1 , 2 , 3 . Stagnation conditions 

at the corresponding stations are designated by subscripts 
si , s2 and s3 . 

The following quantities are arbitrarily specified or 

known at stations 1 and/or 2 , namely, 

Molecular weights: and W 2 

Ratios of specific heats: snd Y 2 

Stagnation pressure of primary jet: P , 

s u. 

stagnation temperature of primary jet: T^^^ 

Static pressures: = P 2 (= Pg) 

Static temperature of primary jet: T^^ (= saturation tempera- 

ture of primary fluid at pressure P^^) 

Static temperature of secondary jet: T 2 (= Tg) 

Mass flow ratio: m^/m^ = m^/(m^ + m 2 ) = x 

Velocity ratio: ^ 1/^2 ~ ^ 

We treat the fluids at stations 1, 2, 3, 4 and 5 as per- 
fect gases with constant specific heats. 

The velocity ratio of the ejector may be developed as 
follows : 



1 < 
IH* 

il 

11 


^1 


Ml 






^sl 


^1 1 


^2 ^ 


^2 


M 2 


1^2 


R 2 


T2 . 


■ “2 
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We can now use t±iis result and other standard perfect 
gas relations to complete the calculation of key thermody- 
namic properties at stations 1 and 2. It is convenient to 
arrange these calculations sequentially as follows ; 



M 



1 




( 2 . 2 ) 



M., = 



's2 



31 



Yl ^1 ^si 






V, ni Y2 R2 ^2 



^2 ■ ^ 2 

1 + 



s2 



_s2) 

’2 J 

P 



Tj (T32/T2) 



s2 



= P 



s2) 



2 P. 



(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 



Before we can determine the actual conditions at station 
3, the outlet of the ejector, it is first necessary to analyze 
two other hypothetical cases as indicated schematically in 
Fig. 2.1. In all three of these cases the flow is taken as 
steady and adiabatic. 

The first case, that shown in Fig. 2.1(a) , involves a 
hypothetical device which receives two gas streams at stations 
1 and 2 and discharges the resulting mixture at station x . In- 
let conditions at 1 and 2 are identical to those of the actual 
ejector. 
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Fig. 2.1 COMPARISON OF ACTUAL EJECTOR WITH IDEAL 
EJECTOR AND WITH HYPOTHETICAL MINIMUM 
LOSS DEVICE. 
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Recall that the adiabatic mixing of two different gases is 
an inherently irreversible process that always involves a cor- 
responding entropy increase. We stipulate that the first case 
is such that the only irreversibility which occurs is that asso- 
ciated with this mixing. We term this mass mixing to distin- 
guish it from another type of mixing considered below which we 
term momentum mixing. 

The second case, that shown in Fig. 2.1(b) , involves an 
ideal ejector which receives two gas streams at stations 1 and 
2 and discharges the resulting mixture at station y. Again, 
the inlet conditions at 1 and 2 are identical to those of the 
actual ejector. This case satisfies, among other relations, the 
idealized one-dimensional momentum equation for frictionless, 
constant area flow. It can be shown that in this case there is 
an overall entropy increase which includes not only the previ- 
ously mentioned effect of mass mixing but also a further increase 
associated with momentum mixing. 

The third case, that shown in Fig. 2.1(c) , represents the 
actual ejector itself which receives input streams at stations 
1 and 2 which discharges the resulting mixture at station 3. 

All three of the above cases refer to constant area pro- 
cesses in the sense that 

Aj^ = Ay = A3 = (A^ + A^) (2.8) 

Once the hypothetical operating conditions at stations x 
and y have been found, it then becomes possible to express the 
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corresponding actual conditions at station 3. This is accom- 
plished through the use of an ejector effectiveness . 

This parameter is defined more explicitly in the later analysis. 
Values of must be estimated by reference to test data on 

ejectors (i.e.. Ref. 6, 7, 8). 



16 



3. Adiabatic Mixing of Two Streams 

The process in the ejector is treated as the adiabatic 
mixing of two streams. The following relations can be shown 
to apply. 



R, 



R- 



R, 



X 





Y 



X 



T 



sx 



A 



X 



8315 Joule/kg°K (3.1) 



R 

Wi 




(3.2) 


R 

W2 




(3.3) 


R = R-, = 

y 3 


X R^ + (1 - X) 


(3.4) 


Yl \ 
(Yi - 1) 




(3.5) 


^2 ^2 
(Y2 - 1) 




(3.6) 


c = c , 

py p3 


= X C , + (1 - x) C „ 
pi p2 


(3.7) 


^y = ^3 = 


C -5 

P3 


(3.8) 


<=p3 - ’' 3 ' 


T = T , 

sy s3 


= Cp3 '== Si Si + '1 ^ %2 S2' 


(3.9) 


A = A, = 
Y 3 


(A 3 + Aj) 


(3.10) 



Notice that the quantities R , C , y , T and A all 

P s 

have the same values at stations x, y and 3. 

This fact may be used to simplify many of the subsequent 
expressions. It is always permissible to substitute subscript 
3 for subscripts x or y on any of these quantities. 
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4. Mass Flow 



It is convenient temporarily to omit station subscripts 
1, 2, 3, X, y and thereby develop certain needed mass flow 
relations in generalized form. Thus the mass flow across an 
arbitrary station may be written and developed as follows. 



fh = 




= PAM 




P .A 
s 



T /T 
s 



(p^/p) 



M 



(4.1) 



This may be rewritten as 



P A 

m = — ^ f(M) 






(4.2) 



where the auxiliary function f(M) is defined as follows 



/t /T 

f(M) - H 



= M 



1 + 






(Y + 1) 
2(y - 1)' 



(4.3) 



Eq. (4.2) can be applied specifically to stations 1, 2, 
X, y, 3 as follows 



P T A. 

m, = ^ - f(M.) 

^ /r t /y ^ 
1 sl'^'l 



(4.4) 



P - A 

m- = ^ — — f(M ) 

/R T 

2 s2' 2 



(4.5) 
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(4.6) 



also 



m = 

X 



^sx ^3 
'^^3 '^53/^3 



f(M ) 

X 



m = 

y 



P A, 
sy 3 



f(M^) 



(4.7) 



m 3 = 



^3 ^ 

V ^3 '^33/^3 



f(M3) 



(4.8) 




= X 



(4.9) 



ih = A =A-=Ai+A- 
X y 3 1 2 



(4.10) 



From Eqs. (4.4), (4.5), (4.9) and (4.10) we readily find 



that 



\ 2 / 

Then from Eq. (2.10) we infer also that 

" TTT'a^TaP' 




^ ""si ^2 ^<»2> 

"^2 '"s2 W 



(4.11) 



(4.12) 



Since all quantities on the right side of Eq. (4.11) are 
now known, Eqs. (4.11) and (4.12) fix the area ratios 
and . 

From Eqs. (4.5), (4.6), (4.9), (4.10) and (2.8) we deduce 
further that 
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f (M^) 



(4.13) 



(1 - A^A ) 
' a - X) 




**2 '^ s 2 ’'^3 



All quantities on the right side of Eq. (4.13) except 

P are now known. The method of determining P is ex- 
sx ^ sx 

plained in a later section; refer to Eq. (5.15). Once P 

s ^ 

has been found, Eq . (4.13) fixes f(M ) whereupon Eq . (4.3) 

fixes itself. The solution of Eq. (4.3) for M when 

f(M) is known involves an iterative procedure which is ex- 
plained later in this section. 

From Eqs. (4.6), (4.7), (4.8) and (4.10) we also find that 



P f(M ) = P f(M ) = p _ f(M-) 
sx X sy y s3 3 



(4.14) 



Assuming P and M known, this relation along with Eq. (4.3) 
fixes P when M is specified; it also fixes f(M_) and 
when P ^2 is specified. The methods of calculating 
and P ^2 ^.re explained later; refer to Eqs. (6.8) and (7.14). 

Fig. 4.1 is a rough sketch of the function f(M) as defined 
by Eq. (4.3). This fianction vanishes at M = 0 and at M = « . 
The maximum value occurs at M = 1 and equals 



max 



2 


1 — 1 
<N 


Y + 1 





(4.15) 



It is clear from the figure that in the range 0<f< t 

Eq. (4.3) has two positive real roots, one si±»sonic and the 
other supersonic. Usually the supersonic root must be discarded 
because it is not consistent with a further constraint imposed 
by the second law of thermodynamics. 
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Fig. 4.1 APPROXIMATE SKETCH OF FUNCTION f(M) 
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In order to solve Eq. (4.3) for M when values of f 

and Y are specified, we utilize an iterative procedure based 

on Newton's method. Let M denote the nth trial value of 

n 

M . We then define a function F(M ) and its derivative 

n 

I 

F (Mj^) as follows. 



F' (M ) 
n' 



F(M ) 
n 



= 



= f 



1 + 



Y - 1 
2 



M 



n 



(Y + 1) 

2Ty"-' 1) 



- M 



n 



(4.16) 





/ \ 


(" 2 ") 


1 + ( ) M ^ 


\ 2 / n J 



(3 - Y) 
2(y - 1) 



- 1 (4.17) 



The (n + l)St trial value of M may now be taken as 

(4.18) 



F(M ) 

^(n -M) = ^n - 



F (Mj^) 



The cycle of calculations defined by Eqs. (4.16), (4.17) 
and (4.18) is repeated until the result for converges to 
a stable value at the desired level of accuracy. 

Thus 



M = ilim M 
n 

n " 



(4.19) 



It is readily apparent that the value of M so obtained 
satisfies Eq. (4.3) for the prescribed values of f and y . 
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5 . Entropy 

Consider the case shown in Fig. 2.1(a). The specific 
entropy s of the discharged gas may be expressed in the 
form 




where and P ^2 the partial pressures of the two com- 
ponents of the gas mixture. Symbols P^ and denote the 

pressure and temperature of the ambient atmosphere. The en- 
tropy of each component is assigned the value zero at this 

reference state P , T . 

o o 

Incidentally, the primary fluid, being condensible, might 
well exist only in the liquid state at the reference condition 
P^ , T^ . Nevertheless, so long as it is in the gaseous state 
at station x , Eq. (5.1) may still be used. 

The partial pressures of the two components are propor- 
tional to the respective mol fractions. Hence 




X 

X + (1 - x)R2 






(1 - x) R 2 
X Rj^ + (1 - x) R 2 




(1 - x)R2 



R 



3 




(5.2) 



(5.3) 
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Upon substituting Eqs. (5.2) and (5.3) into (5.1) and 



simplifying we may obtain the result in the form 




whe re 



(5.4) 



^3 ~ ^ (x Rj^) - ( 1 - x)R 2 X^n [(1 - x)R 2 ] 

(5.5) 



It can be shown that this quantity s^^ represents the 
entropy increase associated with the mass mixing of the two 
different gases. 

Moreover, since the entropy of the stagnation state s 
is by definition identical to that of the corresponding static 
state X , we may replace T and P in Eq. (5.4) by T 
and P , respectively. Thus we obtain 



s = 

X 



m 



+ C 



p3 





(5.6) 



In previous studies of this series, we have considered 
only the special case in which the same fluid is employed for 
both the primary and secondary streams. In that case we have 



R 



1 




(5.7) 



C,=C-=C, = C 
pi p2 p3 p 



(5.8) 



^1 = ^2 = ^3 = Y 



(5.9) 



Moreover, in this special case there is no mass mixing, so 
that we must set 
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0 



(5.10) 



s 

m 

Consequently, Eq. (5.6) now simplifies to 

= 0 + Cp in - R in (5.11) 

Next reverting to the case of two different gases, let 

At denote a , time interval during which unit mass crosses 
station x, . During this same interval x units of mass cross 

station 1 and (1 - x) uinits of mass cross station 2. The 

entropy which leaves across station x must equal the entropy 
which enters across stations 1 and 2 plus the entropy increase 
caused by the mass mixing. Thus 



s = s + xSt+(1- x)s~ 
X m 1 2 



(5.12) 



where 



Si = Zn 



®2 ^p2 



"A 




In 1 


'!£1 


^ 0 1 


1 




! si \ 


- 


£n 




/ 


2 


\'’o 



(5.13) 



(5.14) 



Upon eliminating s between Eqs. (5.6) and (5.12), we 

find that s also cancels from the result. It is then a 
m 

simple matter to solve for the exit stagnation pressure in the 
form 



(^) = k ^p3 (^) - =1 - (1 - ==> = 2 ! = ' 
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Thus the sequence defined by Eqs. (5.13), (5.14) and 

(5.15) now fixes P . Next referring back to Eq. (4.13), 

we can calculate the value of f (M ) ; finally, from Eqs. 

(4.16) , (4.17), (4.18) we can calculate itself. Once 

P , T == T , and M are known, it is a simple matter to 

calculate the corresponding conditions P , T from the 
usual perfect gas relations. 
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6 . Momentum 



Consider the idealized one dimensional constant area 
ejector shown in Fig. 2.1(b). The momentxam equation for this 
system may be written 

(Pi ^ Pi ''l^l^i * (Pj * P2 ' (Py Py 

The mass conservation relation for this system can be 
expressed in the form 



Pi 



(1 - x) *^y y ^y 



( 6 . 2 ) 



Let us divide the three terms of Eq. (6.1) by the corres- 
ponding three terms of Eq. (6.2) . Notice that the areas can- 
cel out of the result. In this way we obtain 



X 



(?! + Pi V^^) 



+ (1 - X) 



(P2 + P2 V2^) 

P2 ^2 



(P -t- p V ^) 

y y y 

P V 

y y 



(6.3) 



Temporarily dropping the station subscript, we next deve- 
lop the typical term of Eq. (6.3) in generalized format as 
follows . 



(P + pV^) 
PV 



[P H- (yrt)m^] _ 

1= /yW m ' V y m /t /t 

R-l S 




(1 + yM^) 






(6.4) 



Observe that the pressure P cancels from the result. 
It is convenient to define the auxiliary function 
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g(M) = — -- A + (3L_^)m^ (6.5) 

(1 + yM^) ^ 

With this notation Eq. (6.3) may be rewritten in the 

form 

X /\ '^sl (1 - X) / ^2 '^s2 ^ 1 7^3 

g(M^) / “glMp” / Y2 g(My) / Y3 ^ ^ 

Rearranging gives 




Since all quantities on the right are known, Eq. (6.7) 

fixes g(M^) . Then itself can be found by inverting 

Eq. (6.5). Fortunately, an explicit solution is possible in 

2 

this case as Eq. (6.5) can be reduced to a quadratic in M 
The result is 



„2 ^ (1 - ± A - 2(y -v ,5 3, 

1 - Yd - 2Yg ) 

The general character of the function g(M) is sketched 

in Fig. 6.1. Notice that g(M) vanishes at M = 0 , reaches 

its peak value g at M = 1 , and decreases toward the 

^max 

limit g as . It can be shown that 

00 



g 

^max 



1 

/2(y + 1) 



g 



00 



1 / (V -TT 

Y /. 2 



(6.9) 



( 6 . 10 ) 
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Fig. 6.1 APPROXIMATE SKETCH OF FUNCTION g(M) 
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Fig. 6.1 discloses the existence of a supersonic root over 
the range ^oo^^'^^max * root should not be accepted until 

a check calculation is made to determine whether the result is 
consistent with the second law of thermodynamics. The subsonic 
root is found by retaining only the negative sign before the 
radical in Eq. (6.8). 

Once has been found from Eq. (6.8) , the corresponding 

value of may be foiond from Eqs. (4.3) and (4.14). This 

suffices to fix all properties at states y and sy . 
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7. Availability 

The entropies and S 2 of the two input streams 

were defined earlier in Eqs. (5.13) and (5.14). The corres- 
ponding steady flow availability fxinctions with respect to an 
ambient atmosphere at pressure and temperature T^ may 

be written 





(T , - T ) - T S, 

sl 0 0 1 


(7.1) 


♦2 = ^2 


(T _ - T ) - T s, 
S2 0 0 1 


(7.2) 



Consequently the total available energy entering the sys- 
tem becomes 



ti<m = X 'I'l + (1 - 

The respective entropies of the streams leaving at stations 
X , y and 3 are 



s 



X 



s 



s 



y 



3 



®m ^3 



m p J 



S + C -5 
m p3 




R3 .n 
R 3 in 
K 3 in 



(7.4) 



(7.5) 



(7.6) 



where s^ has been previously defined in Eq. (5.5). 

The corresponding availabilities of the streams leaving 
at stations x, y and 3 may therefore be written 
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(7.7) 



ip = C,(T,-T)-T s 
p3 ' s3 o ox 



i/ = C,(T,-T)-T s 
p3 s3 o o y 



= C , (T , - T ) - T s- 
^3 p3 s3 o o 3 



(7. 8) 
(7.9) 



Now consider the losses of availability listed below. 
By substituting the above expressions for the quantities on 
the left and simplifying we readily obtain the expressions 



('P^ - 



X 



('J^^ - 



X 



Thus 








(ri) - 


' ^ ) 


= T s 
O 


m 


- ip ) 

y 


II 

U) 


T 2,n 1 
o ’ 


V ^sy 


- ^^ 3 ^ 


= ^3 


T iln 
o 





(7.10) 



(7.11) 



(7.12) 



These results are very significant. Eq. (7.10) defines 
the loss of available energy caused by the mass mixing of the 
two different gases. Eq. (7.11) defines the further loss of 
available energy caused by the momentum mixing in an idealized 
ejector. Eq. (7.12) defines the augmented loss of available 
energy caused by the momentum mixing in the real ejector. 

It is useful to postulate an empirical relation between 
these last two losses which is shown below on the left. The 
expression on the right then follows from Eqs . (7.11) and 

(7.12) . Thus 



( 4 ^ 

Tf 



X 



X 




= 






(7.13) 
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We term the ejector effectiveness. Its value, 

which is always less than unity, must be estimated from appro- 
priate test data on ejectors. 

If we treat as known, Eq. (7.13) fixes • "^he 

solution is simply 

) 

We can next find f(M 2 ) from Eq. (4.14) and from 

Eqs. (4.16), (4.17) and (4.18). This suffices to fix all pro- 
perties at states 3 and s3 thereby couple ting „ the analysis of 
the ejector. 
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8 . 



Diffuser 



It is convenient to analyze the diffuser in terms of the 

concept of availability in steady flow. However, we wish to 

deal with the characteristics of the diffuser itself and these 

cannot depend on the arbitrary values which happen 

to characterize the condition of the ambient atmosphere. We 

can achieve our objective by defining the availability with 

respect to a reference state / "^3 which characterizes 

the diffuser rather than with respect to the ambient state 

P , T . 
o ' o 

Under these circumstances we may write the appropriate 
availabilities at diffuser inlet and outlet as follows 



These equations maJce use of the fact that 
Also, states s3 and 3 are at the same entropy so that 




( 8 . 1 ) 




( 8 . 2 ) 




(8.3) 



From Eqs. (8.1) and (8.3) we conclude that 




(8.4) 
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This result shows that the available energy at the 
diffuser inlet is simply the inlet kinetic energy itself. 

The loss of availability through the diffuser may be 
found by subtracting Eq. (8.2) from (8.1). The result is 



*''^s4 " '^s3^ = AiJ; = + T^ Jin 



3 3 



s 3 
s4 



(8.5) 



We now define the diffuser effectiveness as 






s4 

S3 



= n, 



( 8 . 6 ) 



so that 



'*^3 - *S4 



_ 



4> 



s 3 



'P 



s3 



(1 - n^) 



Rj Tj in 

Cp3 <^s3 - ^3> 



(8.7) 



This reduces readily to 



Jin 



(1 - n^) 



’s4j 



^3 M3* 



( 8 . 8 ) 



Solving for the pressure ratio gives 



s4 

’s3 



= e 



- (1 - n^) 



Y3 M3^ 



(8.9) 



This is the result required. It fixes when P ^3 , 

M 3 and rij 3 are specified. 
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Assiaming the area ratio A^/A^ to be arbitrarily spec- 
ified, we can find from the continuity relation. Refer- 

ring to Eq. (4.8) we may write 

m, = ^ — - — f(Mj = ^ ^ f(Mj = fti, (8.10) 

/«3 ^s3/V3 ' 

Consequently 

' te) ft) "‘“3’ 

Eq. (8.11) fixes f(M^) • Then follows in the usual 

way from Eqs. (4.16), (4.17) and (4.18). The result suffices 
to fix all properties at states 4 and s4 . Thus 





1 + 




( 8 . 12 ) 



= T V(T -,/Tj 
s3' s3^ 4 



(8.13) 



^3 




= Ps4/(^s4/^4) 



(8.15) 



The configurations studied in earlier work did not incor- 
porate a diffuser. For the purpose of comparing the results 
of the present analysis with corresponding earlier work, it is 
desirable to be able to eliminate the effects of the diffuser 
in certain cases. This can be accomplished in the present 
analysis simply by setting 
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1 



(8.16) 



n 



D 



and 




= 1 



(8.17) 



It is easy to confirm that under these circiams tances 



^s4 


= ^s3 


(8.18) 


^s4 


= ^s3 


(8.19) 


M4 


= M3 


( 8 . 20 ) 


^4 


= A 3 etc. 


( 8 . 21 ). 



and all diffuser effects disappear. 
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9. 



Electrical Power Section 



The power conversion process in an EHD duct is usually 
treated on the basis either of constant area or of constant 
static state. Inasmuch as electrical power output per unit 
mass is small, the numerical differences between the results 
computed by these two methods is negligible. For definite- 
ness in this analysis, however, we assume constant area. 

The electrical power that can be obtained from an EHD 
duct of constant area, negligible change of density and optimum 
length can be estimated from the one dimensional version of 
Poisson's equation which governs the electrical field. The 
solution is well known and will not be derived here; a detailed 
derivation may be found in Ref. (3) . The essential result may 
be written in the form 

where 

Pg = gross electric power output, watts 
e = permittivity of medium 

= 8.854 X 10 farad/m (for any gas) 

Ej^ = dielectric strength of medium, volts/m 
= area of duct (constant) 

= velocity 
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According to the test data analyzed in Ref. 4, the’ dielec- 



tric strength is well approximated by the expression 




(9.2) 



where is a characteristic constant of the medium. The 

data show that Eq. (9.2) applies to air or steam up to about 
lO atmospheres pressure with 



In this study we assxame that Eq. (9.2) can be extrapolated 
up to about 100 atmospheres pressure. We also assume that an 
expression of the same form applies to other media besides air 
or steam but that each medi;am has its own characteristic value 
of the breakdown constant C„. 

In the two-fluid system, the separate values of the 
breakdown constants Cg^ and Cg 2 are usually known, but the 
breakdown constant Cg^ of the resulting mixture is seldom 
known. In the absence of adequate test data bearing on this 
point, we tentatively ass\ame that Cg^ can be estimated from 
the hypothetical relation 



Upon substituting Eq. (9.2) into (9.1) and dividing 
through by the mass flow rate, we find the electrical work out- 
put per unit mass of fluid in the form: 



Cg 2 - 9.49 X 10^ m^ “K/cmb 



(9.3) 




(9.4) 
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By use of standard perfect gas relations, Eq. (9.5) can 



be changed to a more useful form. In this connection set 



(T^4 - T^5) = AT^ 



and 



B = 



^ S3 S 



(9.6) 



(9.7) 



where 



Pq = ambient pressure, N/m 
= ambient temperature, °K 



Notice that S represents a dimensionless version of 
the electrical breakdown constant. 

Eq. (9.5) now yields the important result 



1 




The corresponding gross electrical work output per unit 
mass of fluid is then simply 

W* = Cp3 ATg (9.9) 

Unfortunately, the known value of e and the typical 
experimentally measured values of Cg are very small. Con- 
sequently, the dimensionless breakdown constant 6 is ex- 
tremely small compared to unity. Hence the drop in stagnation 
temperature AT^ and the gross electrical work per unit mass 
W*^ are also typically very small. This is a very basic and 
serious limitation on the performance that can be obtained from 
an EHD generator. 

Treating the process 4-^5 as isentropic, we may write 
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(9.10) 



^3 

(Y3 - 1) 



P c = P y. V m 

s5 s4 \T 



•s5 

s4/ 



By analogy with Eq. (4.8) we write the continuity rela- 
tion as 



*4 = ■ f(^4) = f(M5) = K 

J^3 '^s4^'^3 /^3 ’^s5^^3 



(9.11) 



It then follows that 



This fixes ^(^ 5 ) whereupon may be found in the 

usual way from Eqs. (4.16), (4.17) and (4.18). However, con- 
vergence should now be veiry rapid since we may take as a 
first approximation 

w C9.13) 



This solution now suffices to fix all properties at 
states 4 and s4 . 
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10. Condenser/Separator 

The primary fluid which leaves the condenser/separator 
at station 6 is assumed to be compressed liquid at known 
static pressure ~ known static temperature 

T = T . The kinetic energy of the liquid at station 6 is 
assumed to be negligible. 

Temperature Tg must be equal to or greater than the 
ambient temperature T^ in order to satisfy the requirements 
of a heat balance on the condenser/separator. No attempt is 
made in this analysis to formulate this actual heat balance; 
instead Tg is simply treated as a given or known quantity. 
However, it may be of interest later to investigate the 
effects of varying Tg on the overall performance of the EHD 
system. 

The secondary fluid which leaves the condenser/separator 
and enters the ejector at station 2 consists primarily of a 
low molecular weight non condensing gas. It also contains some 
secondary vapor, but the mass fraction of this vapor is treated 
as negligible in this analysis. Thus the secondary fluid at 
station 2 is treated as dry gas. The kinetic energy at station 
2 is not negligible and must be taken into account in the ana- 
lysis . 

In analyzing the thermodynamic process in the condenser/ 
separator, it is again useful to employ the concept of avail- 
ability in steady flow. The availability function is always 
defined with respect to some suitable reference state. Nor- 
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mally this is fixed by the conditions P , T of the ambient 

o o 

atmosphere. In some cases, however, we seek to characterize 
the performance of some component of the system in a manner 
which is independent of the ambient conditions . In 

such cases it is useful to choose as a reference state some 
condition that actually occurs in the component in question. 

We used such a procedure to good advantage in analyzing the 
diffuser and in defining the diffuser effectiveness • We 

now use it again in analyzing the condenser/separator and in 
defining a friction loss coefficient c^ which characterizes 
the operation of this component. For this purpose we choose 
as the appropriate reference state of zero availability the 
static conditions Pg , Tg in the condenser/separator. 

It also clarifies matters further if we choose the con- 
dition P^ , T- as the datum state of zero enthalpy and zero 
6 o 

entropy. In this way we ensure that our analysis in this sec- 
tion makes no reference whatever to the state P^ , T 

o o 

Remember also that Pg = P 2 = and Tg = . 

Thus we may now write 



S- = s + C - Zn 
5 m p3 



■s5' 



- R^ Jin 



s5 



( 10 . 1 ) 



= 



C ~ Jin 
P2 



’s2 



= 0 



- R, 



s2 



= 0 



( 10 . 2 ) 

(10.3) 



The quantity s^ in Eq. (10.1) is the same as previously 
defined in Eq. (5.5). It expresses the entropy increase asso- 
ciated with mass mixing. 
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The availabilities at stations 5, 2 and 6 are therefore 



^5 = S3 


- 




(10.4) 


to 

II 

0 

to 


(^s5 - 


T 2 ) - 0 


(10.5) 


*6 ‘ ^ ■ 


S 


= 0 


(10.6) 



The irreversibility I of the process which occurs in the 
condenser/separator with respect to the specified reference 
state Pg , Tg may be expressed in terms of the above avail- 
abilities as follows 

I = 4'5 - (1 - x) (10.7) 



We postulate that this irreversibility can also be corre- 
lated with the inlet kinetic energy in terms of an empirical 
friction loss coefficient c^ as follows 



V 



I = c. 



/Cf Y3 R3 



M 5 



( 10 . 8 ) 



where as usual 




(10.9) 



Of course the above friction loss coefficient c^ must 
be estimated from suitable experimental data. 

Assuming that c^ is a specified constant, the trial 
solution obtained for an arbitrary trial value of x and 
arbitrary fixed values of the various other input parameters 
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will not necessarily satisfy Eq. (10.7). It is therefore 

convenient to represent the residual error in this equation 

in the form below. The factor T- is inserted on the left 

3 6 

to make E(x) dimensionless. Thus 

R 3 Tg E(x) = ^5 - (1 - x) ’i '2 - 0 - I (10.10) 



The foregoing relations can be assembled and rearranged 
as follows. Let 




Then the final energy balance equation reduces to 



E(x) = K + K, X (10.13) 

o 1 

By holding all other input parameters fixed and allowing 
X to take on a range of possible values, it will normally be 
possible to find by trial a value of x such that 

E(x) = 0 (10.14) 

This then is the value of x which is consistent with 
the other specified input parameters and which defines a pos- 
sible EHD thermodynamic cycle. 
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Of course, if the various fixed input parameters are not 

properly chosen, Eq. (10.14) cannot be satisfied for any trial 

value of X which lies in the useful range 0<x<l and no 

solution exists under these conditions. 

If a solution exists it can be found by some suitable 

iteration technique. One useful method is to set E(x) = 0 

in Eq. (10.13) and then solve this equation for x but to 

relabel this as x . Then 

new 



K 

X = - 

new 



(10.15) 



Assuming that conditions have been so specified that a 
solution for x actually exists, and assuming that an initial 
trial value of x has been chosen that is sufficiently close 
to the true root, repeated application of Eq. (10.15) should 
finally yield convergence to the true root itself. 

Incidentally, any trial solution for which E(x)<0 re- 
presents a physical impossibility because it violates the 
second law of thermodynamics. Any trial solution for which 
E(x)>0 satisfies the second law but represents a condition 
in which the losses in the condenser/separator exceed the 
value prescribed by Eq. (10.8). Only the case for which 
E(x) = 0 satisfies both Eq. (10.8) and the second law. 
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11. Cycle Efficiency 

Once a solution has been found such that E(x) = 0 , it 
becomes a s traightfoward matter to calculate the correspond- 
ing overall cycle efficiency. 

Because of the pressure drop through the boiler/super- 
heater , the pressure at boiler inlet is slightly 

higher than the stagnation pressure P , of the primary 

s ^ 

fluid at the ejector inlet. We assume that the r3tio P , /P ^ 

si s7 

is a specified constant. Hence with P^^ specified, P^^ 

is also known. Also the pressures Pg = P^ are known. 

* 



Consider the ideal pump work w done per unit mass of 

P 

primary fluid. For the present consider the hypothetical case 
of a reversible pump. Also note that density changes across 
the pxamp are negligible. Hence we may write 



* 



w 



P 




( 11 . 1 ) 



The ideal gross electrical work output per unit mass of 



mixed fluid has earlier been established in Eq. (9.5). It is 

* 

denoted by symbol w^ . 

Let us now denote the net useful electrical work output 



per Tonit mass of primary fluid by symbol may be 

★ -k 

related to w^ and w^ in the following way 



w 



net 




Joule/kg 



( 11 . 2 ) 



where n / the excitation efficiency, allows for the small 
electrical power expended to excite the system and where ri_ r 
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the pump efficiency, allows for irreversibilities in the pump. 

The heat input in the boiler/superheater per unit mass 
of primary fluid may be written 



w 



q. = (h^, - h.) - 

si 6 rip 



(11.3) 



where hg is the enthalpy of the primary liquid at conditions 
, T^ as listed in suitable tables of properties. 

O 0 

The overall thermal efficiency of the cycle can novv be 
calculated from the simple formula 



n 



c 



w 



net 



‘in 



(11.4) 



Both n and w , are useful parameters which charac- 
c net 

terize the overall thermodynamic performance of the EHD sys- 
tem. 

Eq. (11.2) shows clearly how the ejector serves as a 
kind of cimplifier which increases the electrical work output 
per unit mass of primary fluid. This effect is shown by the 
presence of the parameter x in the denominator of the first 
term. Notice the beneficial effect of a low value of x on 
parameters '^riet '^c ‘ Unfortmately , the value of x 

cannot be stipulated independently in advance; it is fixed 
by the other specified input parameters as explained in the 
previous section. 

* 

Notice that owing to the fact that w^ is normally 
very small, it is possible in some circumstances for w . 

^ net 
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to take on values which are actually negative. When this 
happens it means that electrical power output is insufficient 
to drive the pump. 

The relationships developed in this and the preceding 
sections make it possible to carry out systematic parametric 
studies of various one- fluid and two- fluid EHD cycles, with 
realistic allowances for the various losses that occur. Such 
studies can establish optimum design parameters and perfor- 
mance limits under various circumstances. These results in 
turn can finally permit informed conclusions to be drawn con- 
cerning the ultimate feasbility of this general type of EHD 
power generator. 
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12. Calculation Sequence for Two-Fluid System 



In this section we summarize the various equations pre- 
viously derived in the approximate order in which they would 
be used in the calculation of system performance. Also listed 
are the initial input parameters whose values must be specified 
in order to start the calculation and various further input 
parameters whose values must be specified in order to proceed 
with various successive stages of the calculation. 



Initial Input Data ; = ^2 “ ^6 ' "^1 ' ^sl ' "^sl ' ^^sl ~ ^1^ 

R = 8315 Joule/kg°K 




-1 



( 12 . 1 ) 




( 12 . 2 ) 




(12.3) 




1 



(12.4) 



Further Input Data: W , = W., , T., = , P (T^) 

■' • V u Z Z D Vo 






w 



(12.5) 




Verify that co<<l 



50 




Verify that M2<1 




Further Input Data; x (trial value) 




R 



+ (1 

^2 ^2 
(Y2 - 1) 






(1 - X) Cp2 



( 12 . 6 ) 

(12.7) 

( 12 . 8 ) 

(12.9) 

( 12 . 10 ) 

( 12 .. 11 ) 
( 12 . 12 ) 

(12.13) 

(12.14) 

(12.15) 

(12.16) 
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(12.17) 



= P3 

^3 “ ^ 3 ) 



^s3 = ^ Si Si + (1 - S 2 S2^ 



f (Mj_) 



= M, 



^1 “ ^ 2 
1 + ^ 



(Yi + 1) 
2 (Yi - 1) 



f (M^) 



= M, 



^2 ~ ^ 2 
1 + -^— M 2 



(Yo + 1) 

2T77^:-rr 



X 



'’s2\ Pi ^sl ''2 #'«2' 



^2 ^ ~ \^sl/ V^2 S 2 ^1 



(12.18) 



(12.19) 



( 12 . 20 ) 



( 12 . 21 ) 



Aj_ (A^/A2) 

A^ ^ (1 + Aj_/A2) 



Further Input Data : 




£n (P /P ) 
sx' o 



sx 



= P e 
o 



( 12 . 22 ) 



(12.23) 

(12.24) 

(12.25) 



(12.26) 
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(12.27) 



(1 - A./AJ /P 2 
f (1 - X ) 



^3 *^ s3 ''^1 
R — f ^ 77 - £(^2) 
sx/ V^2 ^s2 ^3 



(Ys + 1) 
2 \ 2 (Y 3 “ 1 > 



max 



J3 



+ 1 



Verify that f(M ) = f<f 

X in^x 

For first approximation set M = M- 
Iterate using Eqs. (12.31), (12.32), (12.33) 



F(M ) = f 
n 






(Y 3 + 1) 

2(Yo - !)■ 



- M 



n 



F’ (M^) = f 



(■^3 M 



M ll-'" X M 

n\ 2 n 



(3 - Y 3 ) 
2(Yo - 1) 



F(M ) 

^(n +1) ^ ^n " F' (Mj^) 



Iteration converges to fix M 



X 



g(M, ) = 



-^-T7 \ 1 ^ ) ""l 



(1 + Yi ) 



M- 



g(M-) = 



(1 + Y2 ^2 ^ 



^2 ~ ^ 2 
1 + (^2 ^ ^2 



( 12.2 8 ) 

(12.29) 

(12.30) 



(12. 31) 



(12. 32) 



(12.33) 



(12.34) 



(12.35) 
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(12. 36) 



g(My) 



= g = 



Ih V3 


. 1 


^Sl + (1 - X) |i 


/^^s2'\ 


1/ ^3 




1 Yi g(M2) V 





M 



2 = 



(1 - 2y 3 g^) - y 1 - 2 (Y3 + 1) 
1 - Y3 (1 - 2Y3 g^) 



(12. 37) 



Eq. (12.37) applies provided that g 7 ^ goo “ ~\ 



(Yo - 1 ) 

— (12.38) 



If 9 - solution reduces to 



(Yo - 1) 



My = 



2 y- 



(12.39) 



Once My is known from Eq. (12.37) or (12.39) 



(Y3 t 1 ) 



f(M ) = M 1 + ^ 

y y 1 2 



(^3 - 1) „ - « 



(12.40) 



P = P 



f(M^) 



sy sx f(My) 



(12.41) 



Further Input Data; n. 



1 /n, 



p o = p 1 
s3 sx \ P 




(12.42) 



f(M3) = f(M^) 



Verify that £(^ 3 ) = f<f 



Max 



(12.43) 

(12.44) 
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(12.45) 



For a first approximation take M = M 

y 

Iterate using Eqs. (12.31); (12.32), (12.33). 
Iteration converges to fix . 







Further Input Data; 



P 



s4 







Verify that f(M^) = 

For a first approximation take 

Iterate using Eqs. (12.31), (12.32) 

Iteration converges to fix M. . 




(12.33) . 




(12.46) 

(12.47) 

(12.48) 

(12.49) 

(12.50) 

(12.51) 

(12.52) 

(12.53) 



(12.54) 
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1 + 



(12.55) 



's4 



M ^ 

2 ^^4 



^4 = Ts4/(Ts4/^4) 



(12.56) 



Y- 



P T ^''^3 

^s4 _ -^s4 

T.1 

4 4 



^4 = Ps4/(^s4/^4) 



(12.57) 

(12.58) 



-12 

Further Input Data ; £ = 8.8 54 x 10 farad/m, , Cg2 / 

3 2 

(Cg = 9.5 X 10 m °K/cmb for air or steam) 



S 3 S Si + (1 “ ^2 S2^ 



(12.59) 



6 = 



^ S3 S 



(12.60) 



AT 



’s4 



= B 



'Y3 - 1 ) 
2Yo . 



s4\/ 



■s4) 



^3 " ^ 2 

1 + ( 2 )S 



-1/(Y3 



- 1 ) 
(12.61) 



, at 

;s5 = (1 - ^) 

• /I s 4 

s4 



(12.62) 



s5 



T , V^3 - 1' 

= ^34 

s4 



(12.63) 
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f(M5) = f = 



(12.64) 




Verify that = f<f 

^ 5 max 

For a first approximation take 

Iterate using Eqs. (12.31), (12.32) 

Iteration converges to fix . 



M, 



(12.33) 



(12.65) 

( 12 . 66 ) 




Further Input Data ; 

s„ = R, S-n R - X R, Jin (x R, ) 
m3 3 1 1 



- (1 - x) R 2 Jin [(1 - x) R 2 ] 




(12.67) 

( 12 . 68 ) 

(12.69) 

(12.70) 

(12.71) 




(12.72) 



(12.73) 
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(12.74) 



Verify that K^<0 



X 



new 



K 

_c 

K, 



(12.75) 



Iterate from Eq. (12.14) until x converges. 

Further Input Data: (Psl^^s7^ ' ^6 ' ^^sl " ^6^ ' '^x ' 






“e = =p3 



(12.76) 

(12.77) 



w 



* = ^^s7 ^6^ 



(12.78) 



★ ★> 
n w w 

W ^ = i-ii-l - -£ 

net \ X rip 



RESULT 



(12.79) 



w 

q. = (h T - h-) - — £ 
^in 'si 6 n 

sr 



.(12.80) 



w 



n„ = 



net 



c q 



in 



RESULT 



(12.81) 
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13. OPTIMUM VELOCITY RATIO 



The calculation sequence summarized in the preceding sec- 
tion reveals that a comple cycle calculation requires us to 
specify numerical values for about thirty input parameters. 

One of these is the velocity ratio (V^/V 2 ) = y . The calcula- 
tion procedure outlined above then permits us to find the cor- 
responding value of the mass flow ratio x , assuming that 
there exists a value of x which satisfies the governing 
thermodynaraic relations for the specified values of the other 
input parameters . 

If we now change the value of y but hold all other in- 
put parameters fixed, we can repeat the above procedure and 
find a corresponding new value of x , if such a value exists. 
It is evident, therefore, that x becomes some definite 
function of y as long as we remain in the domain where a 
real solution exists. Of course all other dependent variables 
of the cycle are also functions of y . In particular, the 
overall cycle efficiency is some definite function of y . 

Moreover, for prescribed values of the other input parameters, 
there will be some definite value of y , let us designate it 
as the optimum value yQp-(- • which yields the greatest value 
of that is possible under the specified input conditions. 

Our basic purpose therefore is to determine for any prescribed 
values of the other input parameters, the value of YQpt 
corresponding value of ’^c^max’ course the values of ^opt 

and other dependent variables are also of interest. 
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Once a ;vorking computer program is available to establish 

values of Y . , / x . and so on for prescribed 

■‘opt c max cpt ^ 

values of the other input parameters, we can then proceed to 
study the effects of changes in these other parameters and to 
search systematically for such values of the most important 
input parameters as will yield the best overall performance 
of the system. 
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14. CALCULATION SEQUENCE FOR ONE-FLUID SYSTEM 



The calculation sequence summarized in Section 12 simpli- 
fies for the special case where the same condensible fluid is 
used for both the primary and secondary streams. This special 
case is simimarized below. 

Initial Input Data ; “ ^2 ~ ^6 ' ' '^2 ~ "^6 

Assuming pressure » and quality are specified, we 

may find from tables of properties of the saturated fluid 

hi = hg - (1 - z^)hjg (14.1) 



Sj^ - Sg - (1 - (14.2) 

Then with s^ fixed and specified, we may also find 

T^^ and h^^ from tables of properties of the superheated 
vapor. 

In addition with Tg specified somewhere in the range 
T^ £ Tg ^ Tj^ , and with Pg = P^^ specified, we may find Hg ^ 
enthalpy of the liquid at state 6 from the tables of properties. 

The primary fluid can now. be approximated by an "equivalent 
perfect gas" whose properties are as follows 



Y 



C 



P 




Jin(Tg^/T^)'^ 

i^n(Psi/Pi) j 



("si 



(T 



si 




-1 



(14.3) 

(14.4) 



R = 




T 




(14.5) 
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The secondary fluid is represented by the same "equivalent 

perfect gas" and is therefore assigned the same values of 

Y t C and R . 

P 

Further Input Data; (V /V^) = y 
' X ^ 



Ml 



M, 



(y - 1 ) 



sT 

uh , 



- 1 



•s 2 l 



(14.6) 



ill m' 




W 2 / ll^2 


(14.7) 


that M 2 < 1 


(14.8) 


*• (^)"i 


(14.9) 


T 2 (T^ 2 /T 2 > 


(14.10) 



‘ 2\ Y/(Y - 1) 

^s2 = P2(W^2) 

Further Input Data; x (trial value) 



= xT, + (l-x)To 
s3 si s 2 



f(M^) = 



f(M 2 ) = 



1 + (^S-^)m 2 



^-(Y + 1)/2 (y - 1 ) 



,-(Y + 1)/2(y - 1) 



1 + (34-X)m2 



(14.11) 

(14.12) 

(14.13) 

(14.14) 

(14.15) 



fe) 



/^sl ^“'2> 



sJV T 32 f(Mi) (1 - X) 



(14.16) 
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(14.17) 



A^l 



A1/A2 

(1 + A^/A2) 



Further Input Data ; (ambient conditions) 



= T^) ■ 



(14.18) 



= Cpln(g)- Mn(^ 



(14.19) 



An 



sx ’ 



1 

R 



rrp 



Cp in 



s3 



X - (1 - x) S 2 



sx 



P e 
o 



iln(P /P ) 
sx' o' 



(14.20) 



(14.21) 



f(M ) 

X 



'max 



(1 /^s2' 




f(M2) 



Verify that f(^x^ = f < 



For a first approximation set ~ M 2 and iterate 

Iterate using Sqs. (14.26), (14.27), (14.28). 



(14.22) 

(14.23) 

(14.24) 

(14.25) 



F(M ) 
n 



n 



L_ M 



= f 1 + ( 



Y 



_ ^ 2\(Y + 1)/2(Y - 1) 



2 n 



- M 



n 



(14.26) 



Y + 1 I y - 1 ~ 1) 

F'(M^) = f( ^ 2 + ( ^ ' -, - -)M. ' 



n 



= M - 



F(M ) 
n 



(n + 1) n F ' (M ) 

n 



-1 (14.27) 

(14.28) 
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Iteration converges to fix M, 



X 



M, 



g(M,) = 



g(M2) 



M 



(1 + yMj_) 
M 

2 

(1 + YM^) 



n/i + 



n/i + 



g(M^) = g = 



X \ rsl (1 - x)] r^s2 



gTMirVT33 g(M2> UT33 



^ (1 - 2Yg ) - n/i - 2(y + 1) 

1 - Y(1 - 2Yg^) 



-1 



(14.29) 



(14.30) 



(14.31) 



(14.32) 



It can also be shown that in tlie special case where 



g = 9 o 



= 1 yi 



- 2 — — / the above solution degenerates to 



Once is found from Eq. (14.32) or (14.33) 



f (M- ) = M 

y y 



sy 



= P 



1 + O - o - -^ )M^ 
2 y 



f (H^) 



- (Y + 1I/2(y - 1) 



sx f(M ) 

y 



(14.33) 

(14.34) 

(14.35) 



Further Input Data; n-. 



. = p 

s3 sx 1 P 

V sx i 



1/n;. 



(14.36) 



f(M3) 



sx 

’s3 



f(M^) 



(14.37) 
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Verify tnar £(^ 13 ) 



f < f. 



(14.38) 



For a first approximation, take 1 (14.39) 

Iterate using Eqs . (14.26) , (14.27) , (14.28) . 

Iteration converges to fix • 








3 J 



(14.40) 

(14.41) 



pf = <^s3/^3> 



Y/(Y - 1) 



^3 ' ^3/<^s3/^3> 



Further Inpub Data; n,-^ / 



^s4 “ ^S3 



- (1 - np) (Y 143/3) 



^«‘3> 



Verify tliat f(Mj) = f < f 



nax 



A. 



For a first approximation take ^ ~| 

Iterate using Eqs. (14.26), (14.27), (14.28) 

Iteration coverges to fix . 

T ■ = T o 

s4 s3 



T, 



1 Y 



(14.42) 

(14.43) 

(14.44) 

(14.45) 

(14.46) 

(14.47) 

(14.48) 

(14.49) 
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= T^4/<T=4/^4> 



(14.50) 



s4 



= (T34/T4) 



Y/(y - 1) 



= Ps4/(^s4/^4) 



(14.51) 

(14.52) 



Note; If it be desired to omit the diffuser, merely set 



Hq = 1 and A^/A^ = 1 . The result i: 



’s4 



^ 3 / ^4 
— rp 

■^s3 

= 



(14.53) 

(14.54) 



= P 



s3 



= P. 



M 4 = M 3 etc 

and all effects of a diffuser disappear from the result. 



(14.55) 

(14.56) 

(14.57) 



Further Input Data; e = 8.354 x 10 



-12 



farad/m for all gases 



Cg= 9.5 X. 10^ m^ °K/cmb for air or steam 



2 

3 O 



(14.58) 



AT 



T 



s4 



so 



= 3 



f y - 1) rs4 



^s4 - 



•34 



1 + 



y - 1 



,-l/(y-l) 






(14.59) 



(14.60) 



so 



^s4(Ws4) 



y/(y - 1) 



(14.61) 



f (M„) 






(14.62) 
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Verify that f(M^) = f < (14.63) 
For a first approximation take M- ~ . (14.64) 
Iterate using Eqs. (14.26), (14.27), (14.28) . 

Iteration converges to fix . 






1 + M? 

2 3 






sb 

^5 






= '’s5/<^’s5/^5> 



Further Input Data ; 



Note: For a single fluid cycle ~ 



K, 



= In 



^s5\ 



Y 

(Y - 1) 



Jin 



(^) - 






(14.65) 

(14.66) 



(14.67) 

(14.68) 




^5 (14.69) 



K 



1 




E (x) = X 



(14.70) 

(14.71) 



K 

X = - -O 

new (14.72) 

Return to Eq. (1^*13) and iterate \antil x converges and 



E(x) =0 . 
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Further Input Data: (P_,/P_-,) , , 

313/ O 

'’s7 ' 

w* = C AT 
e P s 



(h3l - hg) , 



'^X ' ’^p 



(14.73) 



(14.74) 



VI 



<'’s7 - 



(14.75) 



w 



net 



n w 
X e 

X 



w* 

"p 



RESULT (Joule/kg) (14.76) 



^in = (h^l - ’-6) - Ir 



(14.77) 



w 



net 



^in 



P£SULT 



(14.78)' 
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15. SAJyiPLE CALCULATION FOR A ONE -FLUID SYSTEM 

In this section we illustrate numerically how the calcu- 
lations for a one “fluid system may be carried out for one 
arbitrary but typical operating point. To find an optimum 
operating condition would require the repetition of similar 
calculations for a large number of such trial operating points 
No attempt is made in this section to find such an optimum, r 
We arbitrarily choose to calculate a case which is in 
some respects similar to that identified as case B in the 
appendix of Ref. (1). Thus, we have for both of these cases - 



Medivim 


Steam 


Maximum pressure 


1500 psia 


Quality at state 1 z 


0.95 


Ejector effectiveness ri 


0.9 


Friction factor c^ 


« — 1 

• 

O 


Dielectric strength 


9,4 90 m? °K/cmb 


Diffuser characteristics 


No diffuser 


Excitation efficiency 


1 


Pump efficiency n 


0.9 


Pressure ratio across P ,/P _ 

, . , si s7 

boxler 


1 



Nevertheless, despite these similarities, these two cases 
are not identical. Case B treats the densities at states 1, 2 
X, y, 3, 4, 5 as equal. The present analysis takes the true 
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density changes at these states fully into account. Also it 
is expedient in the present compressible analysis to define 
the parameters and c^ somewhat differently than they 

are defined in Ref (1) . 

Ref • (1) also shows that a comparatively good performance is 
obtained by choosing 

?! = P 2 = Pg = 600 psia (15.1) 



and 




4.25 



(15.2) 



Therefore we arbitrarily choose the same values for the 
present nunierical example, despite the fact that it is not 
entirely equivalent to case B of Ref (1) . We also choose to 
set Tg = (15.3) 

The above conditions now suffice to fix the following 
properties. These are obtained from Keenan and Kayes steam 
tables and are listed on page lOA of Ref. (1) . 

T T = 1115.9 °R 

si 

Tj_ = ^2 = Tg = 946.2 °R (15.3) 



(h^^ - h^) = 78.9 Btu/lbm 

(hg^ - hg) = 673.9 Btu/lbm 

VJe designate ambient conditions as follows 

Pq = 14.696 psia = 1.013 x 10^ N/m^ (15.4) 

Tq = 520 °R = 2 88. '9 
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The properties of the "equivalent perfect gas" may now be 



computed as follows 



-1 



Y = 1 - 



1 - 



£n(Tg3_/T^)’ 

IFTp^Tpp- 



£n(1115.9/946.2) 
£n(1500/600) 



-1 



= 1.2196 



(15.5) 






(^1 - ^l) 



78.9 



(^sl - ^l) 



(1115.9 - 946.2) 
= 0.4649 Btu/lbm 



(15.6) 



R = 



_ (Y - 1) 



Y 



S 1.2196 (0*4649) 



0.2196 



= 0.08371 Btu/lbm “R = 65.15 ft Ibf/lbm °R (15.7) 

The following constants will prove useful for later appli- 
cation 

1 



^ ~ ^ = 0.1098 



= 1.1098 



(Y - 1) 



= 4.5537 



Y 



(y - 1) 



- 5.5537 



(Y 1) _ c n‘i‘^7 

"TTFiT ■ 



i j,— . Jf.l = 4 0537 
2(y-l) 



(15.8) 



Also the dielectric breakdown constant can be expressed in 
dimensionless form as 



6 = 



2 

eC P 
B^o 



(8.854 X lO”^^) (9.49 x 10^) 



(1.013 X lO''"^) 



(288.9) 



= 0.9678(10) 



-3 



(15.9) 
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The Mach numbers at the ejector inlet are 



Ml = 



(Y - 1) 
1.2781 



(^-1) 






0.2196 '94.62 



1115.9 



- 1) 



(15.10) 



M2 = ( 



^2 \ 



2 _ 1 
% 4.25 



(1) (1.2781) 



= 0.3007 



(15.11) 



The following functions are calculated for later use 



f (M^) = Mj_ (1 + 



Y 



- 1 2 ■ 1)/2(y - 1) 



2^1 ) 



=12781(1 + 0.1098(1.2781)^)” 5.50537 



= 0.5552 



(15.12) 



f (M 2 ) = M 2 (1 + - 2 - - - M2^) 



= 0.3007 (1 -t- 0.1098 (0.3007)^) 5.0537 



= 0.2861 



(15.13) 



max 



^Y + 1^ 
= 0.5907 



2 (Y + 1)/2(y - 1) 



5.0537 



= ( 



2.2196 



) 



g(M, ) = 



J 1 * M J 

(1 + Y M^^) 



1.2781 



1 + 0.1098(1.2781) 



(1 + 1.2196(1.2781) ) 



(15.14) 



= 0.4639 



(15.15) 



g(M2) = 



Ml ^ 



1 ^ " 2 ^ 

(1 + yM2^) 
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0.3007 



1/7 



+ 0.1098(0.3007) 



(1 + 1.2196 (0.3007) ^ ) 
0.2722 



(15.16) 



Stagnation conditions at station 2 are calculated next. 



(T^j/Ti) = 






(Ps2/P2> ' 



s2 



1 + 2 ^ ^^2^ = 1 + 0.1098(0.3007)^ 

1.0099 

T2(Ts2/T2> = 946.2(1.0099) 

955.6 

~ = ( 1 . 0099 )^*^^^^ 

1.0562 

Ps(Ps2/^2^ = 600.0 (1.0562) 

633.7 psia 



(15.17) 



(15.18) 



(15.19) 



(15.20) 



In this particular case the entropies at stations 1 and 2 
happen to be equal. They can be calculated either from the stag- 
nation temperatures and pressures or from the static temperatures 
and pressures. The latter choice gives: 



Si = S2 = CpJln(T^/TQ) - R 

= 0.4649 X.n(946. 2/520.0) - 0.08371 Jin (600 . 0/14 . 696) 

=-0.0322 Btu/lbm °R (15.21) 

Further Data Input : x = 0.472 (trial value) 

The above value of x was found by previous trial calcula- 
tions not shown here. The calculations which follow verify that 
this value of x does indeed satisfy the energy balance equation 
of the system to within the round-off error. 
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Ts 3 = X + (1 - x)Tg2 = 0.472(1115.9) + 0.528(955.6) 
= 1031.3 



(1 - X) 'P 



P o 

(^) 



sl 



T^l £(Mj) 



s2 



f (Mj_) 



(15.21) 



0.472 ,633.7 , 

' 1 c n n r\f 



0.528 '1500.0' 



1115.9(0^j = 0.2103 



955.6 '0.5552' 



(15.22) 



A, 



A^/A2 



0.2103 



(l+AjTAp^ 1.2103 



= 0.1738 



(15.23) 



Since = S 2 in this case, Eq. (14.20) simplifies to 



^n(p^) = ^[c £n(= — ) - s,] 

o ^ o 



P = P 
sx oe 



0.08371 
4.1876 



[0.4649 £n ( b - Q| - -;- -) + 0.0322] = 4.1876 (15.24) 



520.0 



= 14.696e = 967.9 psia (15.25) 



(1 - A,/A3) P 2 

= (1 - X)' ■ (p^) 

SX 



*<” 2 ' 



s2 



_ (1 - 0.1726) ,633.7, 

0.528 '967.9^ 






1031, 3 



955.6 



(0.2861) = 0.3045 (15.26 



Solving for M from Eqs . (4.16), (4.17), (4.18) gives 



M = 0.3225 

X 



(15.27) 



For an ideal ejector 



g(My) = 



g(Mj_) UT 



sl (1 - x) 



s3 



g(M^) 




-1 



0.472 



0.4639 



1115.9 ^ 0.528 



1031.3 0.2722 V1031.3 



-1 



955.6 



= 0.3418 (15.2 



M 



(1 - 2yg^) - y 1 - 



2(y + l)g^ 



1 - Y(1 - 2Yg ) 
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_ [1 - 2.4392 (0 .3418)^] - ^1 - 4. 



4392 (0.3418) 

1 - 1.2196 [1 - 2.4392(0.3418)^] 



= 0.1658 



My = 0.4071 



(15.29) 



M 



f(M ) 

y 



V 77 p (Y + d/2(y - 1) 

[1 ^ L 1 



M ] 

y 



0.4071 



[1 + 0.1098 (0.1658) ] 



5.0537 0-3716 



(15.30) 



sy 



0 3045 

= P_.. V = 967.9 ( 7^^1 ^ ) = 793.1 psia 



sx f (M ) 

y 



0.3716^ 



(15.31) 



^sv 793 1 ^/O • 9 

'wtr* = 

S 3 



f(M 3 ) = (^)f(Mj^) = (ii^) (0.3045) = 0.3799 
s3 



(15.33) 



Solving for from Eqs. (4.16), (4.17), (4.18) gives 

M 3 = 0.4183 (15.34) 

In the absence of a diffuser we may write 



^s4 


= = 1031.3 °R 

S3 


^s4 


= P ^3 = 775.7 psia 


M 4 


= M 3 = 0.4183 


f(M4) 


= f(M 3 ) = 0.3799 



The drop in stagnation temperature through the working 
section is given by the following expressions. 



(15.35) 
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1 



2 



AT 

(^) 

^s4 



•^s4 



[1 + ( - ~ 2 — ) ^ 4^1 



1/(Y - 1) 



o QC-70 /I n\ ~3 ,0 . 1098^ ,775.7 ^ , 520. 0^ 
0.9678(10) ( i. 2196 > <IT:696> <IimT7> 



(15.36) 

-3 



[1 + 0.1098 (0.4183) 



4.5537 



= 1.0722(10) 



AT. 



-3 



(15.37) 



= Tg4(ATg/Tg4) = 1031.3 (1.0722) (10) = 1.1057 *R 



Notice that this temperature drop is only about 1.1 °R or 
about one tenth of one percent of the absolute stagnation tempera- 
ture I 



T c = T . - AT^ = 1031.3 - 1.1 - 1030..2 °R 
S5 S4 S 



(15.38) 



s5 



= Ps4(^s5/^s4) 



y/(y - 1) 



_ ,1030.2, 

- <T03T3^ 



5.5537 



= 771.1 psia 



(15.39) 



fdXI^) = 



"^s5 ^s4 

(^)f(Mj 

^s4 ^s5 ^ 



= 1 



1030.2 ( 775.7 ^ (0.3799) 



1031.3 '771.1' 



= 0.3820 



(15.40) 



Solving for from Eqs. (4.16), (4.17), (4.18) gives 



Me 



0.4210 

rY-l^M 2 _ 



(15.41) 



(T ./Tj.) = 1 + ( '' „ )Mc = 1 + 0.1098 (0.4210) = 1.0195 (15.42) 

So D Z D 



^5 = = rrml = 



(15.43) 



76 



To verify the energy balance in the condenser/separator we 



calculate 





- 5-5537 te( 



1030.2 

946.2 



0.12196 ,1010.5 
2 ^946.2 



) (0.4210) 



2 



0.2330 



(15.44) 



= 1 (— ^ - 1) = 5 5537 ( 

1 (Y - 1) ^Tg 



T 



1030.2 

946.2 



K 



1) = 0.4930 



(15.45) 



The final energy balance now gives 



E(x) = X 



0.2330 + 0.4930 (0.472) 



0.0003 



(15.46) 



This verifies that the energy balance is satisfied to within 
round-off error. It also confirms that the value x = 0.472 is 
the mass flow ratio that corresponds to the various independent 
input conditions stipulated in this sample calculation. 

Now that the value of x has been confirmed, we may complete 
the thermodynam.ic cycle. 

Since the pressure drop across the boiler is taken as negli- 
gible in this example 

= 1500.0 psia (15.47) 

The gross electrical work output per unit mass is 

w * = C AT = 0.4649(1.1057) = 0.5140 Btu/lbm (15.48) 

e p s 

The ideal pump work is 
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* 



w 



p 







p 



6 



144 (1500.0 - 600.0) 

778.26 49.75 



= 3.3472 Btu/lbm 



(15.49) 



The net useful work of the cycle is therefore 
* * 

^x'^e _ (1) (0.5140) 3.3472 

'^net X n 0.472 " 0.9 

P 

= - 2.6301 Btu/lbm (15.50) 



Notice that the net work of the cycle happens to be negative 
at this particular operating point. In other words the electrical 
power output is insufficient to drive the piimp under these condi- 
tions . 

The net heat input to the cycle is 

* 



<lin = - >'6> - ^ - "3.9 - 



3.3472 

0.9 



■= 670.2 Btu/lbm 



(15.51) 



Finally, the overall thermal efficiency of the system at this 
operating point is 



n = '^net _ -2.6301 

q. 670.2 

^in 



0.003924 



(15.52) 



Although the performance at this arbitrary operating conditions 
is very poor, the fact remains that there exists some (as yet 
unknown) value of the parameter ^ 2 ~ ^ which the efficiency 

reaches its maximum possible value. 



78 



This optimum value of y can be determined by further 
systematic trial calculations similar to the above sample 
calculation. 

The foregoing calculations have established the values of 
Tg , Pg and M at all stations. The corresponding values of 
T , P , p and s at these locations can then be determined 
from the following relations. 



^ = 1 + ~ ' ^ )M^ = 1 + 0.1098M^ 

T = Tg/{Tg/T) “R 



P^ Y/(Y - 1) 5.5537 

— = (— ) = (— ) 

P 'T ^ ^T 



p = Pg/(Pg/P) psia 

" = If = f 

■'s 

s = C Zn (A - Mn(pS-) 



T P 

= 0.4649Jin(^|^) - 0 . 08371£n Btu/lbm °R 



(15.53) 

(15.54) 

(15.55) 

(15.56) 

(15.57) 

(15.58) 

(15.59) 



The calculation of the above properties has been completed 
in this way and the results are summarized in Table 15.1. These 
same states are shown on a temperature entropy diagram in Fig. 15.1. 
Notice from the tabulation that in passing through the 

3 

ejector, the fluid density p changes from 1.402 Ibm/ft to 
1.524 Ibm/ft^ , an increase of about 8.5% . In passing through 

3 

the working section, the density changes from 1.524 Ibm/ft to 
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Fig. 15.1 TEMPERATURE- ENTROPY DIAGRAM OF THERMODYNAMIC STATES 



3 

1.515 Ibm/ft , a decrease of well under 1% . These data support 
the principle that it is permissible to neglect the density changes 
across these components for purposes of simplifying the analysis. 
This is particularly justified because at present the uncertainity 
in the value of the empirical ejector effectiveness Ht, far 
exceeds the small error introduced by neglecting these minor density 
variations. Moreover, the present sample calculation applies to a 
velocity ratio 4.2 5 which is an extreme case. For 

efficient ejector operation the ratio only 

slightly greater than unity and under these circumstances the den- 
sity change across the ejector is expected to be far less than 
8% . Thus there is little substance in the argument that an 
analysis which neglects these minor compressibility effects must 
be seriously in error. Such an argument runs counter to the 
fundamental principles of fluid mechanics. A more correct appraisal 
of the situation is that the inclusion of these detailed compressi- 
bility effects can be expected to improve the accuracy of the 
analysis to some degree, but that they should not be expected to 
change the overall results in any radical way. 
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16. Nomenclature 



A 


“ 


area 


a 


- 


speed of sound 




- 


characteristic breakdown > 




- 


specific heat at constant 




- 


friction loss coefficient 


Eb 


- 


breakdown electric field 


E(x) 


- 


Eqn. 10.13 


F(M^) 




Eqn. 4.16 


f (M) 


- 


Eqn. 4.3 


^max 


- 


Eqn. 4.15 


g(M) 


- 


Eqn. 6.5 


g 


- 


Eqn . 6.7 


g 

^max 


- 


Eqn. 6.9 


^00 


- 


Eqn. 6.10 


h 


- 


enthalpy 


I 


- 


irreversibility 




- 


Eqn. 10.11 


^1 


- 


Eqn. 10.12 


M 


- 


Mach number 


m 


- 


mass flow rate 


^e 


- 


gross electrical power 


P 


- 


pressure 




- 


heat input 


R 


- 


gas constant 


R 


— 


universal gas constant 
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I 



s 



entropy 



s 



m 



T 

V 

W 



* 



w 



e 



* 



w 

w 

X 

X 



p 

net 

new 



z 



entropy of mixing 

temperature 

velocity 

molecular weight 

gross work out per unit mass 

ideal pump work per unit mass 

useful electrical output per unit mass 

ill ^ / m ^ 

Eqn. 10.15 
V;^/V3 

quality of vapor 



B 

Y 

A( ) 
e 



CO 



Eqn. 9.7, dimensionless breakdown constant 
ratio of specific heats 
change of 

permittivity of the medium 

overall cycle efficiency 

Eqn. 7.13, ejector effectiveness 

Eqn. 8.6, diffuser effectiveness 

pump efficiency 

excitation efficiency 

gas density 

steady flow availability 

mass ratio of condensible primary to dry secondary 
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Subscripts 



0 - ambient conditions 

1 - primary, entrance to ejector (exit from boiler/superheater) 

2 - secondary, entrance to ejector (exit from separator) 

3 - mixture, exit from ejector entrance to diffuser 

4 - mixture, exit from diffuser entrance to generator 

5 - mixture, exit from generator entrance to condenser 

6 - primary, exit from condenser entrance to pump 

7 - primary, exit from pump entrance to boiler 

f - saturated liquid 

G - gas or secondary 
g - saturated vapor 
n - next value 
s - stagnation 
V - vapor 

X - exit station of hypothetical minimum loss device 
y - exit station of ideal ejector 
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